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Abstract 

The  physical  behavior  of  refractive  fluid  interfaces  in  turbulent  compressible  separated 
shear  layers  and  the  resulting  structure  of  the  aero-optical  wavefronts  have  been  inves¬ 
tigated  using  new  experimental  studies  and  physical  modeling.  The  experiments  utilize 
the  new  AFOSR/UCI  variable-pressure  facility  which  generates  large  Reynolds  number 
flows  with  elevated  pressures  in  the  test  section  crucial  for  direct  imaging  of  the  refractive 
field  and  interfaces.  Direct,  non-intrusive,  and  non-integrated  imaging  of  the  refractive 
index  field  in  purely  gaseous  flows  is  achieved  using  laser  induced  fluorescence  of  acetone 
vapor  molecularly  premixed  in  air  and  recorded  with  a  high-resolution  intensified  digital 
camera  system,  simultaneously  with  profiles  of  the  propagated  laser  wavefronts  recorded 
with  a  high-resolution  Shack-Hartmann  sensor.  Three  key  accomplishments  have  been 
made.  Firstly,  direct  examination  of  the  aero-optical  interactions  at  every  location  along 
the  laser  propagation  path  in  the  flow  has  been  achieved.  Secondly,  a  new  aero-optical 
interaction  mechanism  has  been  discovered  in  which  non-monotonic  behavior  of  the  cu¬ 
mulative  aero-optical  aberrations  can  occur  depending  on  the  location  of  the  dominant 
refractive  interfaces.  Thirdly,  significant  resolution  robustness  has  been  observed  for  the 
effects  of  the  refractive  field  on  the  laser  aberrations.  These  three  accomplishments  are 
aimed  toward  the  development  of  new  directed  energy  capabilities  for  airborne  vehicles. 
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Catrakis,  Haris  J. 
Aguirre,  Roberto  C. 
Nathman,  Jennifer  C. 
Garcia,  Philip,  J. 
Zubair,  Fazlul  R. 
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Graduated  with  M.S.  in  November  2005,  UC  Irvine 
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Ph.D./M.S.  Graduate  Student,  UC  Irvine 


AFRL  Interactions  and  Industrial  Interactions: 

The  PI  maintains  contact  with  the  following  members  of  the  Directed  Energy  (DE) 
Directorate  at  AFRL,  Kirtland  AFB,  in  Albuquerque,  NM: 


Martinez,  Rudy  Senior  Scientist  &  DE  Strategic  Planner 

Schlic,  Vern  Senior  Scientist,  Laser  Technology 

The  AFRL  interactions  focus  on  basic  high-energy  laser  and  low-energy  laser  aspects 
of  aero-optical  interactions  for  high-speed  and  high- maneuverability  aircraft,  such  as 
tactical  fighter  aircraft,  for  the  development  of  hardware  such  as  novel  laser  turret  designs 
with  flow/beam  control  as  well  as  software  such  as  computational  simulation  codes  and 
airborne  directed-energy  simulators. 

Industrially,  the  PI  maintains  contact  with  Dr.  Michael  Jones  of  Lockheed  Martin  (Ft. 
Worth,  TX)  regarding  airborne  directed-energy  aspects  and  Dr.  James  Trolinger  of 
MetroLaser  (Irvine,  CA)  regarding  remote  laser  imaging  of  turbulence  properties  for 
airborne  directed-energy  systems. 
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2.  Description  of  Major  Research  Accomplishments 


In  aero-optics,  as  discussed  in  review  articles  [1,  2],  one  of  the  basic  goals  is  to  relate 
physically  the  behavior  of  turbulence-generated  refractive  index  fields  to  the  structure  of 
turbulence-distorted  optical  wavefronts.  Ever  since  some  of  the  earliest  studies  in  the  field 
of  aero-optics  [3,  4],  advances  in  both  the  understanding  of  turbulence  and  its  relation  to 
the  optical  behavior  have  been  sought  theoretically,  experimentally,  and  computationally 
(e.g.  [5,  6,  7,  8,  9,  10]).  Knowledge  of  the  aero-optical  interactions,  identification  of 
the  dominant  flow  mechanisms  in  the  generation  of  the  aero-optical  distortions,  and 
optimization  of  the  aero-optical  interactions  by  using  flow  control  where  possible  and 
optical  control,  are  needed  in  both  basic  and  applied  problems  in  aero-optics  especially 
where  laser  beam  propagation  needs  to  be  utilized  [1].  Improved  knowledge  of  laser 
interactions  in  turbulent  flows  is  useful  for  a  variety  of  fundamental  and  technological 
problems  ranging  from  basic  methods  for  probing  the  behavior  of  large  Reynolds  number 
flows  in  the  laboratory,  in  the  field,  and  for  remote  sensing  [11,  12,  13],  to  practical 
techniques  for  beam  propagation  from  directed  energy  systems  on  airborne  vehicles  [1,  7, 
10,  14],  for  example.  At  large  Reynolds  numbers,  the  flow  structure  and  resulting  laser 
beam  behavior  are  accessible  at  present  only  through  experiments  either  in  the  laboratory 
or  in  the  field.  Furthermore,  one  can  appreciate  that  knowledge  of  the  propagated  beam 
behavior  only,  or  of  the  integrated  flow  behavior  only,  does  not  uniquely  relate  to  the 
aero-optical  interactions.  This  is  because  the  optical  aberrations  resulting  from  refractive- 
index  variations  integrate  the  flow  structure  along  the  beam  propagation  path.  Thus  one 
needs  a  means  to  enable  the  examination  of  the  aero-optical  interactions  at  each  location 
along  the  beam  path. 

With  regard  to  refractive  turbulent  fields,  the  basic  physical  interactions  between  the 
propagating  optical  wavefronts  and  the  aberrating  turbulent  flows  occur  across  the  re¬ 
fractive  fluid  interfaces  [10] .  Thus,  an  improved  understanding  of  the  turbulent  interfacial 
behavior  is  needed  [15,  16,  17,  2,  18].  At  large  Reynolds  numbers,  however,  turbulent 
refractive  interfaces  are  highly  irregular  and  span  a  wide  range  of  scales.  Thus,  ex¬ 
perimental  studies  at  large  Reynolds  numbers  are  particularly  desirable  as  a  means  to 
complement  computational  studies  and  physical  modeling.  This  is  because  direct  numer¬ 
ical  simulations  are  at  present  limited  to  moderate  Reynolds  numbers  [19]  and  large-eddy 
simulations  can  be  sensitive  to  sub-grid  scale  modeling  [20,  8].  Furthermore,  depending 
on  the  problem  of  interest,  one  would  like  also  to  incorporate  compressibility  effects  and 
flow  geometry  effects  on  the  modeling  of  the  large  Reynolds  number  behavior.  Previous 
work  in  aero-optics  includes  studies  of  both  the  large  scale  properties  [21,  6,  9,  10]  and 
small  scale  properties  [22,  23,  8].  The  available  studies  indicate  that  the  near  field  aero- 
optical  effects,  in  turbulent  shear  flows,  are  dominated  by  the  large  scale  flow  structures 
(e.g.  [9,  10]).  In  addition,  knowledge  of  the  large  scale  aero-optical  behavior  can  be 
expected  to  be  useful  regarding  modeling  of  the  cumulative  contributions  of  the  small 
scales,  relative  to  the  large  scales,  for  various  Air  Force  applications  involving  directed 
energy  propagation. 

In  the  research  conducted  in  this  grant  period,  a  key  focus  is  on  the  role  of  refractive 
interfaces  in  the  large  scale  aero-optical  interactions  in  high  Reynolds  number  separated 
compressible  shear  layers.  For  this  purpose,  the  new  variable  pressure  flow  facility  at 
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UC  Irvine  developed  by  the  PI  with  prior  AFOSR/DURIP  support  (F49620-02- 1-0277) 
is  utilized  to  generate  the  separated  shear  layers.  The  freestream  consists  of  air  seeded 
with  acetone  vapor  and  the  ambient  gas  is  air.  The  flow  imaging  technique  consists  of 
direct  imaging  of  the  refractive  index  field  and  interfaces  using  laser  induced  fluorescence 
and  digital  photography  techniques.  An  ultraviolet  laser  sheet  excites  the  acetone  vapor 
fluorescence  which  is  in  the  visible  spectrum  and  is  recorded  with  a  digital  camera  system. 
This  enables  recording  of  two-dimensional  spatial  images  of  the  refractive  index  held.  The 
Reynolds  number  based  on  the  visual  thickness  is  Re  ~  6  x  106  and  the  convective  Mach 
number  is  Mc  ~  0.4.  The  refractive  index  fields  are  analyzed  to  identify  the  refractive 
interfaces.  Complex  networks  of  spatially  isolated  high  gradient  boundaries  are  evident 
both  in  the  interior  of  the  how  and  in  the  outer  regions.  Examination  of  the  data  shows 
that  these  high  gradient  boundaries  contain  crucial  information  regarding  the  production 
of  large  scale  aero-optical  aberrations.  The  extent  to  which  the  aero-optical  aberrations 
depend  on  reduction  in  refractive  held  information  is  also  investigated  in  the  context  of 
physical  modeling  based  on  the  high  gradient  interfaces. 

The  aero-optics  variable  pressure  how  facility  utilized  in  the  present  investigations  is 
shown  in  figure  1.  The  facility  permits  experiments  at  elevated  test  section  pressures, 
in  the  range  1  <  p  20  atm,  in  order  to  enable  higher  signal-to-noise  ratio  how  imag¬ 
ing  than  would  be  possible  at  atmospheric  or  sub-atmospheric  pressures.  For  imaging 
techniques  based  on  laser  induced  fluorescence,  such  as  the  method  presently  used,  the 
resulting  image  signal-to-noise  ratio  is  greater  due  to  the  increased  local  density  of  air 
molecules  as  a  result  of  the  higher  operating  pressures. 

The  laboratory  facility  primarily  consists  of  a  blow-down  wind  tunnel  contained  inside 
a  variable  pressure  vessel.  The  main  aero-optics  vessel,  in  which  the  tunnel  is  housed,  is 
oriented  vertically  and  is  shown  on  the  right  part  of  figure  1.  The  gas  reservoir,  which 
supplies  fluid  to  the  tunnel,  is  oriented  horizontally  and  is  shown  in  the  left  part  of  figure  1. 
The  reservoir  is  rated  at  a  maximum  supply  pressure  of  ~  3,  000  psi  or  ~  200  atm.  The 
aero-optics  vessel  that  houses  the  test  section  is  rated  at  a  maximum  operating  pressure 
of  ~  300  psi  or  ~  20  atm.  The  main  vessel  is  eight  feet  in  height  and  four  feet  in  diameter, 
enabling  two  persons  to  stand  upright  in  the  interior  when  reconfiguring  the  tunnel  and 
its  test  section.  The  main  vessel  has  a  quick-release  entrance  door  of  two  feet  in  diameter 
that  facilitates  access  in  and  out  of  the  vessel  in  order  to  reconhgure  as  necessary  the  blow¬ 
down  wind  tunnel  and  its  test  section  contained  in  the  interior  of  the  vessel.  Both  the 
main  pressure  vessel  and  the  reservoir  vessel  are  “U” -stamp  ASME  rated.  Two  ellipsoidal 
cap  sections,  for  each  vessel,  are  welded  at  the  two  ends  of  the  main  cylindrical  sections. 
Both  vessels  including  the  quick-release  door  mechanism  were  fabricated  by  Melco  Steel. 
The  main  vessel  and  the  reservoir  vessel  are  equipped  with  NPT  fittings  for  attachments 
of  vacuum  lines,  pressure  gauges  and  transducers,  and  other  instrumentation  necessary 
for  flow  management. 

Extensive  optical  access  to  the  interior  of  the  main  pressure  vessel  is  available  through  five 
high  optical  quality  spectrosil  windows  each  of  10-in.  diameter  and  3.25-in.  thickness. 
Four  of  the  five  windows  are  vertical  optical  ports  on  the  four  sides  of  the  vessel.  The 
fifth  window  is  a  horizontal  optical  port  on  the  top  of  the  vessel.  Combined  with  laser 
diagnostics,  as  shown  in  figure  2,  the  facility  is  capable  of  quantitative  high  resolution 
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Figure  1:  Photograph  of  the  aero-optics  variable  pressure  flow  facility  at  UC  Irvine 
developed  by  the  PI  with  prior  AFOSR/DURIP  support.  The  main  pressure  vessel  is 
visible  in  the  right  part  of  the  photograph.  It  has  a  large  interior  and  cross  section, 
extensive  optical  access,  and  a  quick  release  door,  to  facilitate  investigations  of  aero- 
optical  phenomena  at  variable  pressures.  Combined  with  laser  and  imaging  diagnostics, 
this  facility  enables  high  resolution  flow  imaging  and  beam  measurements  of  aero-optical 
flows  at  large  Reynolds  numbers  with  compressibility  effects.  The  blow-down  wind  tunnel 
and  test  section  are  housed  inside  the  main  vessel. 
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Figure  2:  Schematic  of  the  laser  illumination  and  flow  imaging  configuration.  A  pnlsed 
ultraviolet  laser  beam  is  shaped  into  a  parallel  laser  sheet,  by  utilizing  a  parabolic  mirror, 
and  is  propagated  through  the  quartz  windows  of  the  pressure  vessel  and  into  the  test 
section,  cf.  figure  1.  A  high  resolution  digital  intensified  charge  coupled  device  (CCD) 
camera  captures  the  visible  fluorescence  from  the  acetone  vapor  that  is  molecularly  seeded 
in  air.  The  laser  sheet  propagation  direction  and  the  flow  imaging  direction  are  mutually 
orthogonal. 
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imaging  of  compressible  turbulent  flows  at  large  Reynolds  numbers. 

The  flow  geometry  presently  examined  is  the  separated  turbulent  shear  layer  which  is  an 
example  of  a  separated  flow  relevant  in  a  wide  range  of  applications  in  aero-optics  (e.g. 
[5]).  A  schematic  of  the  flow  geometry  is  shown  in  figure  3.  In  the  actual  experiments,  the 
flow  is  directed  vertically  upward  and  the  laser  propagation  is  horizontal,  cf.  figure  2.  The 
separated  shear  layer  represents  a  basic  flow  configuration  of  interest  for  future  studies  of 
more  complex  flow  geometries.  For  example,  while  separated  flows  over  curved  surfaces 
have  an  unsteady  separation  point,  the  sharp  90°  corner  of  the  present  flow  geometry 
ensures  a  fixed  separation  point.  Studies  of  the  unforced  separated  shear  layer  can  provide 
useful  baseline  information  to  evaluate  the  extent  to  which  active  flow  control  techniques 
such  as  large  scale  disorganization  or  regularization  can  be  effective  (e.g.  [14,  1]). 

In  contrast  to  two-stream  shear  layers  which  typically  have  smaller  growth  rates  (e.g. 
[24]),  the  separated  shear  layer  or  single-stream  shear  layer  has  a  relatively  large  growth 
rate.  A  large  growth  rate  is  desirable  in  laboratory  studies  because  it  facilitates  higher 
resolution  flow  imaging  due  to  the  large  flow  scales.  While  this  is  preferable  for  fun¬ 
damental  studies  such  as  the  present  ones,  one  needs  to  keep  in  mind  that  in  practical 
aero-optical  systems  a  reduction  in  the  growth  rate  is  generally  desirable  since  that  would 
limit  the  extent  of  the  aberrating  flow  field.  It  is  essential  to  appreciate  that  the  sep¬ 
arated  shear  layer,  as  indicated  in  figure  3,  does  not  exhibit  reattachment  as  long  as 
the  perpendicular  wall  around  the  separation  point  extends  sufficiently  to  ensure  ample 
entrainment  of  the  ambient  fluid  [25].  This  is  in  contrast  to  flows  over  backward-facing 
steps,  of  finite  step  height,  which  exhibit  re-attachment. 

In  the  present  experiments,  the  elevated  test  section  pressure  (p  ~  3  atm)  and  the  rela¬ 
tively  large  size  of  the  vessel  have  resulted  in  a  Reynolds  number  of: 

Re  ~  6  x  106  ,  (1) 

based  on  the  visual  thickness  L  ~  0.1m  of  the  shear  layer  in  the  far  field,  where  the 
laser  induced  fluorescence  imaging  is  conducted,  and  on  the  freestream  velocity  which 
is  U  ~  300  m/sec  for  a  convective  Mach  number  of  Mc  ~  0.4.  The  high  Reynolds 
number  achieved  corresponds  to  flow  conditions  well  above  the  mixing  transition,  i.e. 
the  change  in  mixing  behavior  and  specifically  the  pronounced  increase  in  the  mixed 
fluid  fraction  that  is  directly  associated  with  the  transition  to  fully  developed  turbulence, 
which  requires  at  least  Re  104  [26], 

Figure  4  shows  examples  of  flow  shadowgraphs  at  convective  Mach  numbers  of  Mc  ~  0, 
Mc  ~  0.2,  and  Mc  ~  0.4,  with  increasing  compressibility  corresponding  to  the  images 
shown  from  top  to  bottom  in  the  figure.  The  shadowgraphs  were  recorded  in  the  present 
studies  by  collimating  10-inch  diameter  beams  with  a  parabolic  mirror  and  propagating 
the  beams  in  a  spanwise  direction  normal  to  the  shear  layers.  The  test  section  pressure 
for  the  present  studies  is: 

p  ~  3  atm .  (2) 

The  choice  for  this  test  section  pressure  was  made  because  it  was  found  to  provide 
sufficient  imaging  signal  for  the  refractive  index  field  as  discussed  below  in  §3.  For  a 


LASER  PROPAGATION 
DIRECTION 

Figure  3:  Schematic  showing  the  separated  free  shear  layer  flow  geometry  and  the 
laser  propagation  geometry.  A  convergent  nozzle  is  employed  with  exit  dimensions  of 
~  7.5cmx30cm  to  produce  a  subsonic  freestream  which  separates  past  the  90°  corner 
shown.  The  perpendicular  wall  below  the  separation  point  extends  sufficiently  vertically 
and  spanwise  to  ensure  ample  entrainment  of  the  ambient  fluid.  This  feature,  combined 
with  zero  downstream  pressure  gradient,  ensures  a  separated  free  shear  layer  with  linear 
growth  rate  and  without  reattachment.  The  larger  and  smaller  rectangles  shown  indicate 
the  extents  of  the  fields  of  view  used  for  the  shadowgraphs  and  for  the  laser  induced  flu¬ 
orescence  images,  respectively,  in  the  present  studies.  In  this  schematic,  the  freestream 
direction  is  horizontal,  from  left  to  right,  and  the  laser  propagation  direction  is  vertical, 
from  bottom  to  top.  The  orientation  in  the  experiments,  with  the  vertical  tunnel,  is  90° 
counterclockwise  to  the  orientation  of  this  schematic,  i.e.  the  actual  flow  orientation  is 
vertical  and  the  laser  propagation  direction  is  horizontal,  cf.  figures  1  and  2. 
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Figure  4:  Whole  field  shadowgraphs  of  separated  shear  layers  with  increasing  compress¬ 
ibility  (from  top  to  bottom).  The  top  image  is  for  a  nearly  incompressible  shear  layer, 
the  middle  image  is  for  a  convective  Mach  number  of  Mc  ~  0.2,  and  the  lower  image  is 
for  a  convective  Mach  number  of  Mc  ~  0.4  with  a  Reynolds  number  of  Re  ~  6  x  106 
based  on  the  visual  thickness  at  the  right  edge.  The  images  show  the  aero-optical  in¬ 
tensity  variations  for  10-inch  diameter  beams  propagated  spanwise  normal  to  the  shear 
layers.  The  freestream  fluid  is  a  purely  gaseous  air/acetone  mixture,  above  each  shear 
layer,  and  ambient  air  is  below  each  shear  layer.  Large  scale  organized  vortices  with  sec¬ 
ondary  braid  vortices  are  evident  in  all  images  of  these  single-stream  shear  layers.  The 
flow  organization  and  the  reduction  of  the  growth  rate  with  increasing  compressibility 
are  consistent  with  previous  studies  of  two-stream  shear  layers  [24,  27]. 


10 


convective  Mach  number  of  Mc  ~  0.4,  corresponding  to  the  shadowgraph  shown  in  the 
bottom  in  figure  4,  the  freestream  Mach  number  is: 

~  0.9.  (3) 

The  growth  rate  L'  =  dL/dx  of  the  compressible  separated  shear  layer  at  Mc  ~  0.4  was 
found  to  be: 

L'  dL/dx 

-  rs-/  - 

-h;nc  hZ/jnc  / dx 

where  L[nc  =  dLinc./da;  is  the  growth  rate  for  the  incompressible  separated  shear  layer. 
This  is  in  agreement  with  previous  studies  of  the  effect  of  compressibility  (e.g.  [28,  29, 
30,  31,  32,  33,  34,  35,  27,  36]). 

The  present  study  focuses  on  the  flow  conditions  with  Mc  ~  0.4  which  is  a  moderate 
level  of  compressibility  in  shear  layers  (e.g.  [27,  37]).  In  the  context  of  aero-optics, 
the  moderate  compressibility  conditions  associated  with  Mc  ~  0.4  are  directly  relevant 
to  practical  flight  conditions  for  airborne  directed  energy  systems  (e.g.  [5,  9]).  In  the 
following  sections,  results  are  presented  on  laser  induced  fluorescence  imaging  of  the 
refractive  turbulent  held  (§3),  the  physical  behavior  of  the  refractive  turbulent  interfaces 
(§4),  and  physical  modeling  of  the  aero-optical  interactions  (§5),  for  Mc  ~  0.4. 

In  order  to  conduct  direct  imaging  of  the  refractive  turbulent  held  interfaces  in  purely 
gaseous  how,  the  present  study  utilizes  acetone  vapor  to  seed  the  air  at  a  molecular  level. 
This  provides  three  signihcant  advantages  in  comparison  to  other  approaches.  Firstly,  the 
use  of  acetone  vapor  provides  a  signihcant  laser  induced  fluorescence  signal  at  elevated 
pressures  [38,  39]  in  comparison  to  weaker  signals  such  as  those  arising  from  Rayleigh 
scattering  in  air  (e.g.  [40]).  Secondly,  for  the  present  how  imaging  experiments,  the  ace¬ 
tone  vapor  is  desirable  because  it  ensures  purely  gaseous  fluid  in  the  test  section  and  thus 
it  provides  a  molecular  marker  in  contrast  to  seeding  with  liquid  droplets,  fog  particles,  or 
smoke  particles  for  example  [41],  Thirdly,  and  most  importantly  regarding  beam  steering, 
whereas  approaches  using  high  Rayleigh  scattering  cross-section  gases  such  as  ethylene 
introduce  highly  undesirable  intensity  streaks  (e.g.  [6]),  the  present  approach  results  in 
refractive  index  held  images  that  are  free  of  intensity  streaks  and  thus  eliminating  any 
signihcant  beam  steering  in  the  turbulent  region  of  the  how.  This  is  achieved  by  utilizing 
a  relatively  low  seeding  level  of  acetone  vapor  for  the  fluorescence  imaging.  The  impor¬ 
tance  of  streak-free  imaging  is  that  one  does  not  have  to  resort  to  approximate  methods 
of  attempting  to  remove  the  streaks  as  was  done  in  previous  work  [6] .  Streak-free  imaging 
of  refractive  helds  provides  a  means  to  gauge  whether  or  not  signihcant  beam  steering 
occurs.  In  the  present  experiments,  there  is  no  evidence  of  intensity  streaks  in  the  laser- 
sheet  images.  This  is  because  care  was  taken  to  ensure  that  sufficiently  low  seeding  of 
acetone  vapor  in  air  was  chosen.  If  there  were  signihcant  beam  steering  occurring  in  the 
turbulent  region  of  the  how,  then  there  would  be  signihcant  convergence  and  divergence 
of  the  laser  rays  and  thus  laser  intensity  streaks  (non-uniformities)  with  higher  intensities 
near  beam  convergence  regions  and  lower  intensities  near  beam  divergence  regions.  Clear 
examples  of  signihcant  beam  steering  and  intensity  streaks  are  evident  in  previous  work 
[6].  Thus,  if  there  were  signihcant  beam  steering,  then  one  would  see  intensity  streaks 
generated  in  the  laser  sheet  as  it  propagates  through  the  refractive  held.  The  absence  of 
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Figure  5:  Compressibility  scaling  of  the  growth  rate  of  the  large  scale  visual  thickness 
for  the  present  work  (asterisk  symbol)  and  for  previous  studies  (symbols  as  listed  in 
legend)  [28,  29,  30,  31,  32,  33,  34,  35,  27,  36].  The  growth  rate  for  the  present  primary 
compressibility  condition  of  Mc  ~  0.4  is  shown  and  is  in  accord  with  previous  studies  of 
compressibility  effects. 
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intensity  streaks  in  the  present  laser  sheet  images,  such  as  figure  6  shown  further  below, 
is  evidence  that  there  is  no  significant  beam  steering. 

In  the  present  studies,  acetone  vapor  molecularly  mixed  in  air  is  generated  by  evaporating 
liquid  acetone.  To  prepare  the  purely  gaseous  air/acetone  mixture,  a  stream  of  air  is 
passed  through  a  liquid  bath  of  acetone  using  the  bubbling  method  [39].  The  air/acetone 
mixture  is  supplied  to  the  reservoir  which  is  pressurized  prior  to  each  experimental  run. 
The  main  vessel  is  pressurized  with  filtered  air  at  p  ~  3  atm  prior  to  exhausting  the 
reservoir  gas  into  the  blow-down  tunnel  and  test  section.  The  acetone  vapor  concentration 
in  air  is  regulated  by  controlling  the  mass  flux  of  air  through  the  bubbler.  For  the  present 
experiments,  an  acetone  vapor  concentration  of  ~  2%  by  volume  is  generated  in  order  to 
ensure  minimal  absorption  of  the  laser  power,  as  shown  in  previous  related  experiments 
on  laser  induced  fluorescence  (e.g.  [16]),  and  in  order  to  maintain  conditions  below  the 
lower  flammability  limit  of  acetone.  Fluorescence  of  acetone  vapor  requires  excitation  by 
ultraviolet  illumination  as  discussed  below. 

As  indicated  in  figure  2,  an  ultraviolet  beam  is  generated  from  a  pulsed  NckYAG  laser. 
The  illumination  source  employed  is  a  moderate  energy  laser  system  (Spectra  Physics 
Model  INDI  40-10)  which  produces  ~  70  mJ  per  pulse  at  a  wavelength  of  266  nm.  This 
ultraviolet  wavelength  corresponds  to  frequency  quadrupling  at  the  fourth  harmonic  of 
the  fundamental  infrared  1064  nm  output  of  the  laser.  The  laser  beam  is  shaped  into  a 
laser  sheet  with  a  width  of  ~  15  cm  and  a  thickness  of  ~  150  /jrn.  A  parabolic  mirror 
ensures  that  the  laser  sheet  that  is  produced  is  of  parallel  extent.  The  laser  sheet  is  prop¬ 
agated  through  the  quartz  windows  of  the  pressure  vessel  and  into  the  high  Reynolds 
number  separated  shear  layer  in  the  test  section.  The  laser  sheet  orientation  corresponds 
to  a  streamwise  slice  of  the  shear  layer,  cf.  figure  3.  The  extent  of  the  laser  sheet  that 
illuminates  the  shear  layer  is  ~  15cmxl5cm.  The  laser  sheet  width  of  ~  15  cm  corre¬ 
sponds  to  the  downstream  extent  of  the  imaged  shear  layer.  The  large  scale  transverse 
extent  of  the  shear  layer  at  the  imaging  location  is  L  ~  10  cm.  The  incident  laser  sheet 
excites  the  acetone  vapor  in  the  air  to  generate  visible  (blue)  fluorescence.  The  elevated 
test  section  pressure  of  p  ~  3  atm  results  in  an  increase  in  the  fluorescence  signal,  rela¬ 
tive  to  atmospheric  or  sub-atmospheric  test  section  pressures,  in  agreement  with  previous 
studies  [39]. 

The  laser  induced  fluorescence  imaging  signal  is  recorded  on  a  high  resolution  digital 
intensified  charge-coupled  device  (CCD)  camera  with  ~  1,  000  xl,  000  pixels.  The  camera 
system  is  oriented  along  a  direction  normal  to  the  laser  sheet  and,  therefore,  perpendicular 
to  the  incident  laser  propagation  direction.  The  digital  camera  system  employed  is  an 
intensified  CCD  camera  by  Stanford  Photonics  (Mega-10Z)  with  enhanced  sensitivity  in 
the  blue  spectrum.  This  enables  the  recording  of  the  two-dimensional  spatial  slices  of 
the  acetone  vapor  concentration  field.  The  fluorescence  signal  images  are  normalized  and 
calibrated  pixcl-by-pixcl  using  post-run  images  recorded  at  uniform  concentration  as  well 
as  pre-run  images  that  correspond  to  pure  air  in  the  test  section.  The  normalization  and 
calibration  procedure  is  similar  to  the  method  described  in  a  previous  study  [16].  In 
the  present  experiments,  the  compressible  turbulent  mixing  of  the  freestream  gas  with 
the  ambient  gas  generates  the  refractive  index  field  corresponding  directly  to  the  acetone 
vapor  concentration  field. 
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An  example  of  a  high  resolution  calibrated  and  normalized  image  of  the  refractive  index 
held  of  the  separated  high  Reynolds  number  shear  layer  is  shown  in  figure  6,  for  the  same 
flow  conditions  of  Re  ~6x  106,  Mc  rsj  0.4,  and  p  ~  3  atm.  A  total  of  twenty  uncorrelated 
images  were  recorded  at  the  same  flow  conditions.  As  indicated  in  the  schematic  of 
figure  3,  the  image  shown  corresponds  to  a  two-dimensional  spatial  streamwise  slice. 
The  image  provides  the  instantaneous  refractive  index  held  since  the  laser  pulse  temporal 
duration  is  <  10  nsec  and  thus  faster  than  the  shortest  time  scales  of  the  flow.  The  colors 
shown  in  the  image  correspond  to  variations  in  the  refractive  index.  Pure  ambient  air 
corresponds  to  the  white  region  located  at  the  bottom  of  the  image.  The  freestream 
fluid  is  visible  as  dark  blue  in  the  upper  region  of  the  image.  The  intermediate  colors 
correspond  to  mixed  fluid.  For  these  images,  the  signal-to-noise  ratio  is  estimated  from 
the  measurements  to  be  in  excess  of  ~  100  :  1  which  is  large  enough  to  enable  accurate 
extraction  of  fluid  interfaces  such  as  refractive  index  interfaces  or  high  gradient  interfaces. 
In  flow  imaging  experiments  of  high  Reynolds  number  flows  [42,  6,  37,  12],  there  are 
resolution  limitations  associated  with  the  finite  number  of  camera  pixels  and  with  the 
Unite  laser  sheet  thickness.  The  present  image  data  can  be  viewed  as  coarse-grained  data 
which  are  useful  for  characterizing  the  large  scale  refractive  index  held  behavior. 

As  is  evident  in  the  quantitative  visualization  in  figure  6,  a  relatively  sharp  interface  is 
present  between  the  freestream  fluid  and  the  mixed  fluid  (upper  part  of  image)  as  well  as 
between  the  ambient  fluid  and  the  mixed  fluid  (lower  part  of  image).  This  is  consistent 
with  previous  investigations  that  have  focused  on  the  outer  fluid  interface  between  pure 
fluid  and  mixed  fluid  (e.g.  [43,  6,  44,  16]).  The  sharpness  of  the  outer  interface  is  also 
consistent  with  theoretical  considerations  of  diffusion  (e.g.  [45]).  The  outer  interface 
is  the  dynamic  interface  whose  location  marks  the  extent  of  the  mixed  fluid  region. 
The  relative  sharpness  of  the  outer  interface  depends,  in  general,  on  both  the  Reynolds 
number  and  the  Schmidt  number.  Although  the  Schmidt  number  for  the  present  studies 
is  of  order  unity,  i.e.  corresponds  to  relatively  large  molecular  diffusion  in  comparison  to 
mixing  in  liquids,  the  high  Reynolds  number  provides  the  high  small-scale  strain  rates 
which  limit  the  local  extent  of  molecular  diffusion  in  the  vicinity  of  the  interface  [44], 
The  pronounced  vorticity/scalar  correlations  present  in  fully  developed  turbulent  mixing 
maintain  the  sharpness  of  the  outer  interface  (e.g.  [16]). 

A  scalar  held,  such  as  the  refractive  index  held  n(x,  t),  can  be  physically  represented 
in  terms  of  its  scalar  interfaces.  In  the  present  case,  the  refractive  index  interfaces  are 
diffusive  interfaces  across  which  molecular  diffusion  occurs.  Also,  in  the  context  of  aero- 
optics,  the  interactions  of  the  optical  wavefronts  with  the  fluid  occur  physically  across 
these  interfaces.  Properties  of  turbulence-generated  isosurfaces  of  scalar  helds  g(x,  t),  i.e. 
g(x,  t )  =  const.,  have  been  examined  previously  for  various  how  conditions  in  the  context 
of  mixing  and  other  applications  (e.g.  [15,  46,  16]).  Associated  with  the  refractive  index 
held  n(x,  t),  the  refractive  isosurfaces  are  given  by: 

n(x,  t)  =  const.  (5) 

Although  the  isosurfaces  are  key  how  elements  for  physical  modeling  of  phenomena  sen¬ 
sitive  to  interfacial  behavior,  the  isosurfaces  have  mathematically  zero  thickness  (e.g. 
[15]).  Thus,  it  is  desirable  to  consider  and  examine  physical  thickness  aspects  of  these 
turbulent  interfaces  (e.g.  [47]). 
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Figure  6:  Laser  induced  fluorescence  example  of  a  high  resolution  direct  quantitative 
image  of  the  refractive  index  turbulent  field  in  the  large  Reynolds  number  compressible 
separated  shear  layer,  cf.  the  flow  geometry  schematic  in  figure  3.  The  imaging  of 
the  refractive  field  in  purely  gaseous  flow  is  achieved  using  laser  induced  fluorescence  of 
acetone  vapor  molecularly  mixed  in  air.  The  Reynolds  number  is  Re  ~6x  106  based  on 
the  large  scale  visual  thickness,  the  convective  Mach  number  is  Mc  ~  0.4,  and  the  test 
section  pressure  is  p  ~  3  atm.  The  freestream  is  from  left  to  right  in  the  upper  region. 
The  colors  shown  correspond  to  the  refractive  index  field  and  have  been  added  to  aid 
the  eye  in  discerning  the  refractive  turbulent  interfaces.  Dark  blue  is  for  the  high-speed 
freestream  (top)  and  white  is  for  the  ambient  air  (bottom). 
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Each  refractive  interface  can  be  thought  of  as  consisting  of  fluid  residing  between  the  two 
neighboring  isosurfaces: 

n(x,  t)  +  dn  =  const.,  and,  n(x,  t)  =  const.,  (6) 


i.e.  corresponding  to  the  refractive  indices  n  +  dn  and  n,  respectively. 


The  local  interfacial  fluid  thickness  can  be  defined  in  terms  of  the  refractive  index  gradient 
magnitude.  From  geometrical  considerations,  the  refractive  index  gradient  Vn  is  normal 
to  the  interface  at  every  point  on  the  interface.  As  long  as  the  magnitude  |Vn|  of  the 
gradient  is  finite,  we  can  define  the  interfacial  thickness  hn  per  unit  n  as: 


Mx,t) 


1 

]Vnf 


dn 

ds 


(7) 


where  s  denotes  distance  along  the  local  normal  to  the  interface.  The  actual  spacing 
between  the  isosurfaces,  cf.  equation  6,  is  therefore  hndn.  This  corresponds  to  the  full 
interfacial  thickness. 


By  general  geometrical  considerations,  we  can  express  the  differential  volume  d V(n)  of 
fluid  residing  between  these  two  neighboring  isosurfaces  as: 
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(8) 


where  A{n)  is  the  interfacial  surface  area,  Kef  is  a  reference  volume,  and  hn  dn  is  the 
mean  spacing  between  the  neighboring  isosurfaces  which  can  be  interpreted  as  the  mean 
interfacial  thickness. 


It  can  be  noted  that,  since  the  values  of  the  refractive  index  field  are  dimensionless,  the 
quantity  hn  has  the  units  of  length.  For  simplicity,  therefore,  we  will  refer  to  hn  as  the 
interfacial  thickness.  As  a  result  of  the  strong  intermittency  of  turbulent  flows  (e.g.  [48], 
[47]),  especially  at  high  Reynolds  number  flow  conditions,  the  local  interfacial  thickness 
can  be  expected  to  be  highly  variable. 

The  interfacial  fluid  thickness  hn  for  the  present  data  was  computed  on  the  basis  of 
equation  7.  Since,  in  the  present  measurements,  only  two-dimensional  spatial  information 
is  available  for  the  refractive  index  field,  only  the  in-plane  component  of  the  full  interfacial 
thickness  can  be  evaluated.  This  is  shown  in  figure  7.  Evaluation  of  the  full  value  of 
the  interfacial  thickness  would  require  the  three-dimensional  refractive  index  gradient 
(e.g.  [49]).  Nevertheless,  the  two-dimensional  information  is  useful  in  the  context  of 
aero-optical  interactions  in  terms  of  characterizing  the  in-plane  optical  phase  distortions 
as  demonstrated  in  previous  studies  (e.g.  [6,  10]). 

The  behavior  in  figure  7  indicates  the  presence  of  high  gradient  interfaces,  denoted  as 
dark  regions  in  the  figure.  These  high  gradient  interfaces  are  highly  irregular  and  are 
present  at  the  outer  boundaries  between  mixed  fluid  and  pure  fluid  as  well  as  in  the 
interior.  These  networks  of  relatively  isolated  high  gradient  interfaces,  in  these  two- 
dimensional  images,  are  consistent  with  previous  observations  of  sheet-like  structures  in 


16 


Figure  7:  Example  of  the  in-plane  components  of  the  interfacial  fluid  thickness  hn  corre¬ 
sponding  to  the  refractive  index  image  shown  in  figure  6.  The  interfacial  thickness  field 
shown  was  computed  based  on  equation  7  for  the  two  in-plane  components.  The  magni¬ 
tude  of  the  interfacial  thickness  is  labeled  with  grey  values.  Darker  regions  denote  higher 
gradient  locally  thinner  interfaces.  Highly  irregular  networks  of  high  gradient  interfaces 
are  evident  both  in  the  interior  of  the  flow  and  at  the  outer  boundaries. 
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three  dimensions  in  the  scalar  dissipation  rate  field  (e.g.  [49,  50]).  It  can  be  appreciated, 
in  general,  that  the  wide  range  of  interfacial  scales  at  high  Reynolds  numbers  prevents 
direct  examination  of  the  entire  range  of  thickness  scales.  At  finite  resolution,  such  as 
with  the  present  data,  the  available  refractive  index  field  corresponds  to  a  coarse-grained 
field  with  information  retained  only  at  or  above  the  resolution  scale  (e.g.  [12]). 

Knowledge  of  the  refractive  index  field  is  useful  both  for  evaluating  the  optical  wave- 
front  distortions  and  for  identifying  the  dominant  fluid  interfaces  in  the  aero-optical 
interactions  (e.g.  [6,  9]).  Evaluation  of  the  wavefront  distortions,  given  information  on 
the  refractive  interfaces,  depends  on  the  flow  regime  and  optical  regime.  For  example, 
near-field  interactions  are  well  described  by  geometrical  optics  while  far-field  interactions 
require  wave  descriptions  based  on  Maxwell’s  equations  (e.g.  [7,  51]).  Also,  for  high 
energy  beam  propagation  there  is  two-way  coupling  between  the  fluid  and  the  optical 
waves,  whereas  for  low  energy  beams  there  is  only  one-way  coupling  between  the  flow 
structure  and  the  beam  distortion  (e.g.  [52,  53,  6]). 

For  the  present  studies,  the  measured  refractive  index  fields  are  utilized  with  the  eikonal 
equation  of  geometrical  optics  to  propagate  computationally  the  optical  wavefronts  in 
the  near  field.  As  long  as  the  smallest  turbulent  scales  are  greater  than  the  optical 
wavelength  scales  (e.g.  [54]),  the  geometrical  optics  description  is  sufficient  for  near-field 
aero-optical  interactions.  In  terms  of  the  refractive  index  field  n  and  the  optical  path 
length  (OPL),  the  eikonal  equation  is: 

|  V  (OPL)|  =  n ,  (9) 

(e.g.  [55]).  Solutions  of  this  equation  along  the  individual  ray  paths  correspond  to  the 
OPL  integral  in  terms  of  the  refractive  index  field: 

OPL  (x,  t)  =  /  n(£,  t)di.  (10) 

J  ray 

The  OPL  integral  corresponds  to  inverting  the  eikonal  equation  9.  The  physical  distance 
along  the  beam  propagation  path  for  each  light  ray  is  denoted  as  t. 

The  optical  wavefronts  can  be  represented  as  isosurfaces  of  the  OPL,  i.e., 


OPL  (x,  t)  =  const. 


(11) 


The  OPL  also  corresponds  to  the  optical  phase  or  wavefront  propagation  time  (e.g.  [1]). 
A  useful  way  to  quantify  the  wavefront  distortions  is  in  terms  of  the  optical  path  difference 
(OPD)  given  by: 

OPD  (x,  t)  =  OPL  (x,  t)  -  OPL,  (12) 

where  OPL  corresponds  to  the  ensemble  averaged  OPL  value.  A  useful  measure  of  the 
wavefront  distortion  is  the  root-mean-square  (rms)  value  of  the  OPD  variations,  i.e.: 


OPD 


rms 


^  |OPD(x,t) 


(13) 
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Figure  8:  Flow-derived  optical  wavefront  phase  profile  evaluated  at  the  upper  boundary 
of  the  field  of  view  of  the  refractive  index  field  image  shown  in  figure  6.  The  streamwise 
distance  is  normalized  in  terms  of  the  large  scale  visual  thickness  of  the  flow.  The  en¬ 
semble  averaged  root-mean-square  optical  path  difference  (OPDrms)  was  found  to  be  in 
agreement  with  independent  Shack-Hartmann  measurements  of  the  aero-optical  aberra¬ 
tions  at  the  same  flow  conditions. 


where  ( )ens.  denotes  ensemble  averaging.  The  ensemble  averaged  OPD  value  for  wave- 
fronts  corresponding  to  the  present  database  of  images  was  found  to  be: 

OPDrms  ~  0.2  /irn  .  (14) 

Independent  Shack-Hartmann  measurements  of  the  aero-optical  aberrations  were  con¬ 
ducted  and  found  to  also  result  in  an  ensemble  averaged  value  of  OPDrms  «  0.2  /irn,  i.e. 
in  agreement  with  the  flow-derived  value  of  OPDrms.  It  must  be  noted  that  this  OPDrms 
value  is  for  the  present  flow  conditions.  In  general,  it  can  be  expected  to  be  a  function 
of  Reynolds  number,  Mach  number,  flow  geometry,  and  other  parameters.  The  OPDrms 
value  can  be  used  to  normalize  the  OPD  profile  for  each  distorted  wavefront.  Also,  this 
quantity  is  directly  related  to  a  measure  of  beam  attenuation  known  as  the  Strchl  ratio 
which  is  a  key  quantity  for  evaluating  the  effects  of  the  refractive  index  field  on  the  peak 
intensity  of  the  propagated  beam  [1,  6]. 

Using  the  present  measurements  of  the  refractive  index  field  n  and  equation  10,  the 
evolution  of  the  wavefront  OPD  through  the  cross-section  of  the  shear  layer  can  be 
examined.  Figure  8  shows  an  optical  wavefront  phase  profile  corresponding  to  the  upper 
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Figure  9:  Results  of  the  evolution  (top)  of  the  optical  path  difference  (OPD)  for  wave- 
fronts  propagating  through  the  refractive  field  (bottom)  of  the  high  Reynolds  number 
compressible  separated  shear  layer,  using  the  image  data  in  figure  6  reproduced  here  for 
comparison.  Successive  OPD  profiles  are  shown  at  several  transverse  distances  through 
the  flow,  with  the  wavefronts  propagating  vertically  from  bottom  to  top. 
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Figure  10:  Expanded  visualization  of  the  top  part  of  figure  9  for  the  cumulative  behavior 
of  the  OPD  for  optical  wavefronts  propagating  through  the  compressible  separated  shear 
layer.  As  stated  for  figure  9,  successive  OPD  profiles  are  shown  at  several  transverse 
distances  through  the  flow,  with  the  wavefronts  propagating  vertically  from  bottom  to 
top.  The  streamwise  distance  is  normalized  by  the  large  scale  visual  thickness  of  the 
flow.  The  uppermost  OPD  profile  was  evaluated  at  the  upper  boundary  of  the  held 
of  view  of  the  refractive-index  held.  The  ensemble-averaged  root-mean-square  optical 
path  difference  (OPDrms)  was  found  to  be  consistent  with  independent  Shack- Hartmann 
measurements  of  the  aero-optical  aberrations  at  the  same  how  conditions. 
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boundary  of  the  field  of  view  of  the  refractive  index  held  image  shown  in  figure  6.  Figure  9 
(top)  shows  successive  profiles  of  the  OPD  at  several  transverse  locations  through  the 
shear  layer,  for  planar  incident  wavefronts.  The  corresponding  refractive  index  held 
is  shown  on  the  bottom  in  figure  9,  reproduced  from  figure  6,  to  facilitate  the  visual 
inspection  of  the  aero-optical  interactions.  The  wavefront  is  propagated  vertically,  from 
bottom  to  top,  in  figure  9  (upper  part).  The  hrst  OPD  prohle,  near  the  bottom  of 
the  sequence,  has  planar  structure  since  it  has  not  yet  entered  the  shear  layer  and  thus 
corresponds  purely  to  the  incident  wavefront.  The  second  OPD  prohle  from  the  bottom 
of  the  sequence  shows  the  initial  distortion  of  the  rightmost  part  of  the  wavefront  as 
it  passes  through  the  sharp  outer  interface  of  the  shear  layer.  Additional  large  scale 
features  imposed  on  the  wavefront  structure  by  the  shear  layer  interfaces  are  evident  as 
the  wavefront  propagates  further  into  the  shear  layer.  These  are  evident,  for  example,  in 
the  upper  three  prohles  in  the  sequence.  Successive  nearly-uniform  large  scale  regions, 
or  plateaus,  in  the  refractive  index  held  separated  by  irregular  high  gradient  interfaces 
appear  to  form  a  primary  mechanism  in  the  observed  aero-optical  interactions.  The 
entire  evolution  of  the  wavefront  shape  is  shown  in  figure  11.  The  observed  behavior 
in  the  OPL  prohle  sequence  demonstrates  the  utility  of  refractive  index  held  imaging  to 
examine  the  role  of  the  refractive  interfaces  in  the  aero-optical  interactions. 


The  structure  of  the  refractive  index  held  and  its  effect  on  the  propagation  of  the  wave- 
fronts,  as  depicted  for  example  in  the  sequence  in  figure  9,  suggests  that  the  dominant 
contributions  to  the  wavefront  distortions  can  be  identified  in  terms  of  the  high  gradi¬ 
ent  refractive  interfaces.  To  do  so,  one  can  utilize  the  physical  thickness  hn(x,  t)  of  the 
interfaces  as  defined  above  in  equation  7. 


The  role  of  the  thickness  of  the  refractive  interfaces  in  the  aero-optical  interactions  can 
be  seen  directly  by  expressing  the  OPL  in  equation  10  as  an  integral  with  respect  to  the 
refractive  index  n  rather  than  the  propagation  distance  £,  i.e. 

OPL(x,  t)  =  f  n(£,t)  hU}£  \dn\ ,  (15) 

-Lay 

where  hnj  is  the  interfacial  thickness  component  in  the  direction  of  optical  propagation: 


h 


n,l 


1 

|Vra|*  ‘ 


(16) 


This  takes  into  account  the  orientation  of  the  interface  relative  to  the  light  ray  direction, 
with  |V/r|£  denoting  the  component  of  the  magnitude  of  the  refractive  index  gradient  in 
the  i  direction,  i.e.  in  the  direction  of  the  optical  ray  propagation,  i.e. 


|Vra|*  = 


dn 

d£ 


(V) 


The  absolute  values  in  equations  15  and  17  are  needed  since  the  refractive  index  n  could 
be  locally  increasing  or  decreasing  as  the  light  rays  propagate. 

Whenever  the  interfaces  are  not  normal  to  the  optical  rays,  the  gradient  magnitude 
component  will  be  smaller  than  the  full  gradient  magnitude  and  therefore  the  interfacial 
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thickness  through  which  the  ray  propagates  will  be  effectively  larger.  If  a  part  of  a 
refractive  interface  is  locally  normal  to  an  optical  ray,  the  gradient  component  |Vn|^ 
will  be  the  same  as  the  full  gradient  magnitude  |Vn|.  Thus,  in  general,  the  gradient 
magnitude  component  must  always  be  less  than  or  equal  to  the  full  gradient,  i.e. 

|Vn|^  =  |Vn|  |cosd|  <  |Vra| ,  (18) 

which  means  that  the  effective  interfacial  thickness  must  always  be  greater  than  or  equal 
to  the  full  thickness,  i.e. 

h  a ./  |  j-  =  hn  |  sec  6*|  ^  hn  .  (19) 

|Vn|  i 

The  angle  6,  taken  as  —  n  <  6  <  tt,  quantifies  the  orientation  of  the  refractive  index 
gradient,  which  is  always  normal  to  the  fluid  interface,  relative  to  the  local  optical  prop¬ 
agation  direction.  It  can  be  noted  that  the  above  relations,  as  written,  require  that  the 
local  gradient  magnitude  is  finite,  i.e.  |Vn|  %  0.  They  can  be  modified,  in  a  straightfor¬ 
ward  manner,  in  any  regions  of  uniform  refractive  index,  i.e.  regions  of  |Vn|  =  0.  For 
example,  hnj  dri  can  be  interpreted  as  a  Dirac  delta  function  whose  integral  is  the  extent 
Ad  of  the  uniform  index  region  in  the  optical  propagation  direction,  and  its  contribution 
to  the  OPL  integral  becomes  A  (OPL)  =  n  At',  in  that  case. 

The  role  of  the  interfacial  thickness  enables  a  physical  model  for  the  large  scale  aero- 
optical  interactions,  by  appreciating  that  the  high  gradient  locally-thin  interfaces  provide 
the  dominant  contributions  to  the  wavefront  distortions.  Figure  12  shows  the  comparison 
between  the  full  OPD  profile  and  a  modeled  OPD  profile  where  a  gradient  threshold  is 
employed  to  retain  only  ~  50%  of  the  original  refractive  index  field  values  that  correspond 
to  local  gradients  above  the  threshold.  For  the  ensemble  averaged  data,  there  is  less  than 
~  10%  difference  in  the  root-mean-square  OPD  at  the  present  flow  conditions.  Because 
the  high  gradient  interfaces  are  spatially  isolated,  as  is  evident  in  the  image  in  figure  7, 
the  low  gradient  regions  in  between  the  high  gradient  interfaces  are  modeled  to  a  first  ap¬ 
proximation  as  uniform  index  regions  using  the  gradient  at  each  boundary  through  which 
each  ray  propagates.  It  is  helpful  to  appreciate  that  although  the  high  gradient  interfaces 
have  a  relatively  small  component  hnj  of  the  thickness,  they  contribute  significantly  to 
the  OPL  integral  in  equation  15  because  the  differential  dn  of  the  refractive  index  is 
relatively  large  across  such  interfaces.  In  this  study,  only  the  in-plane  components  of  the 
interfacial  thickness  are  available  and  thus  figure  12  corresponds  to  the  in-plane  optical 
wavefront  distortions  only.  The  out-of-plane  component  would  be  needed  to  capture  the 
out-of-plane  wavefront  distortions. 

The  present  approach  enables  an  examination  of  the  effect  of  filtering  the  refractive  field, 
i.e.  by  retaining  only  high  gradient  interfaces,  on  the  optical  wavefront  phase  aberrations 
for  a  wide  range  of  filter  values.  Figure  13  shows  the  manner  in  which  the  difference  in 
the  root-mean-square  OPD,  between  the  full  wavefront  phase  profiles  and  the  modeled 
wavefront  phase  profiles,  depends  on  the  amount  of  reduction  in  the  refractive  field  infor¬ 
mation  as  a  result  of  the  filtering.  The  percentage  values  of  refractive  field  information 
reduction  correspond  to  the  extent  of  the  filtered  regions  of  the  refractive  field,  i.e.  the 
regions  excluded  as  a  result  of  modeling  by  retaining  only  high  gradient  interfaces.  Higher 
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Figure  12:  Modeled  OPD  profile  (dotted  curve)  compared  to  the  full  OPD  profile  (solid 
curve),  corresponding  to  the  wavefront  propagated  through  the  separated  shear  layer 
in  figure  9.  The  modeled  OPD  profile  is  computed  using  a  modeled  refractive  index 
held  consisting  of  high  gradient  interfaces,  with  low  gradient  regions  replaced  by  uniform 
index  regions,  where  the  gradient  threshold  employed  retains  only  ~  50%  of  the  original 
refractive  index  held  values.  For  the  ensemble  averaged  data,  a  difference  of  ~  10%  is 
observed  in  the  root-mean-sqnare  OPD,  between  the  modeled  wavefront  profiles  and  the 
full  wavefront  profiles,  at  the  present  flow  conditions. 
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Figure  13:  Dependence  of  the  difference  in  the  root-mean-square  OPD,  between  the 
full  wavefront  phase  profiles  and  the  modeled  wavefront  phase  profiles,  on  the  amount  of 
reduction  in  the  refractive  field  information.  The  percentage  of  refractive  field  information 
reduction  corresponds  to  the  extent  of  the  regions  of  the  refractive  held  modeled  by 
retaining  only  high  gradient  interfaces.  Higher  percentages  of  reduction  in  refractive 
held  information  correspond  to  higher  gradient  thresholds  below  which  the  refractive 
held  is  filtered.  The  difference  of  ~  10%  in  the  root-mean-square  OPD  for  a  gradient 
threshold  resulting  in  ~  50%  reduction  in  refractive  held  information  corresponds  to  the 
prohles  shown  in  hgnre  12. 
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percentage  values  of  reduction  in  refractive  field  information  correspond  to  higher  values 
of  the  gradient  threshold  below  which  the  refractive  field  is  filtered.  Significant  robust¬ 
ness  is  evident  whereby  less  than  ~  5%  difference  in  the  root-mean-square  OPD  results 
for  reductions  of  up  to  ~  40%  in  the  refractive  field  information.  For  reduction  percent¬ 
age  values  of  up  to  ~  50%,  less  than  ~  10%  difference  in  the  root-mean-square  OPD  is 
evident.  This  robustness  can  be  understood  by  appreciating  that  it  is  the  high  gradient 
interfaces  that  produce  the  dominant  aero-optical  interactions.  The  observed  robustness 
provides  a  key  ingredient  that  can  be  utilized  in  physical  descriptions,  modeling,  as  well 
as  optimization  studies  in  aero-optics. 

To  gain  further  insight  into  the  one-to-one  aero-optical  correspondence  as  a  function  of 
distance  z  along  the  beam  propagation  path,  we  consider  the  cumulative  optical  path 
length.  In  other  words,  we  consider  the  cumulative  integral  in  terms  of  the  refractive- 
index  field,  as  a  function  of  the  beam  propagation  path  distance: 

OPL(z,f)  =  f  n(£,t)d£  =  f  n(£,t)  hnje\dn\ ,  (20) 

J  path  J  path 

where  hn/  denotes  the  interfacial  thickness  component  in  the  ray  direction  along  the 
beam  path. 

We  can  now  consider  a  cumulative  Strchl  ratio  in  terms  of  the  cumulative  OPDrms,  i.e. 
in  terms  of  the  cumulative  root-mean-square  of  the  optical  path  difference  as  follows: 

SR(-M)  =  Ibif,  (21) 


so  that: 


SR(z,  t)  ~  exp 
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This  cumulative  Strehl  ratio  is  a  useful  quantity  that  can  enable  examination  of  the  one- 
to-one  aero-optical  correspondence  and  cumulative  effects  along  the  beam  propagation 
path. 


Figure  14  shows  examples  of  results  on  the  cumulative  aero-optical  behavior,  i.e.  as  a 
function  of  distance  along  the  dynamic  beam  propagation  path.  The  plot  in  the  top  part 
of  figure  14  shows  the  cumulative  OPDrms  as  a  function  of  propagation  distance  for  an 
individual  flow  realization.  This  is  a  non-monotonic  function  with  increasing  propagation 
distance  for  the  instantaneous  flow  realization.  The  behavior  reflects,  physically,  the 
result  of  different  parts  of  the  wavefront  accumulating  distortions  at  different  rates.  The 
plot  in  the  bottom  part  of  figure  14  depicts  the  cumulative  Strchl  ratio  as  a  function  of 
propagation  distance.  The  non-monotonicity  again  reflects  the  result  of  different  parts 
of  the  wavefront  dynamically  accumulating  distortions  that  can  collectively  decrease  as 
well  as  increase  the  Strehl  ratio,  cf.  equation  22.  The  observed  behavior  in  figure  14  is 
a  new  key  feature  that,  as  can  be  appreciated,  can  only  be  detected  by  examining  the 
cumulative  aero-optical  interactions  along  the  beam  propagation  path.  In  other  words, 
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Figure  14:  Top:  Cumulative  r.m.s.  optical  path  difference  OPDrms,  as  a  function  of 
propagation  distance,  for  dynamic  laser  sheet  propagation  through  a  large-Reynolds- 
number  separated  shear  layer.  Bottom:  Cumulative  Strchl  ratio  SR  as  a  function  of 
propagation  distance  corresponding  to  the  r.m.s.  optical  path  difference  shown  at  the 
top. 


knowledge  of  the  integrated  behavior  only  would  have  been  insufficient  to  discover  this. 
It  is  the  direct,  i.e.  non-integrated,  flow  imaging  that  enabled  the  present  finding. 

The  present  observation  and  finding  provide  a  key  new  aspect  of  aero-optical  interactions 
for  the  aperture-integrated  beam  behavior,  i.e.  for  finite-aperture  beams.  The  extent  to 
which  this  occurs  in  other  flow  realizations  and  at  other  flow  conditions  is  an  important 
future  direction  of  research.  The  study  of  this  will  be  facilitated  with  the  interfacial-fluid- 
thickness  (IFT)  framework  pioneered  by  the  PI  in  previous  AFOSR  work  [10]  to  identify 
the  role  of  the  high-gradient  locally-thin  interfaces  in  generating  the  dominant  contri¬ 
butions  to  the  wavefront  distortions.  The  present  results  provide  new  insight  into  the 
mechanism  by  which  one-to-one  aero-optical  interactions  along  the  entire  beam  propaga¬ 
tion  path  produce,  cumulatively,  the  aberrations  of  the  propagated  laser  beam.  This  new 
mechanism  provides  a  key  ingredient  which  can  be  expected  to  facilitate  the  development 
of  aero-optical  optimization  techniques,  including  active  forcing  by  plasma  excitation  to 
modify  the  primary  refractive  interfaces  in  order  to  reduce  directly  the  beam  aberrations. 

In  the  research  reported  in  this  grant  period,  an  additional  key  focus  is  on  multiresolution 
aspects  of  the  aero-optical  interactions.  The  motivation  is  twofold.  Firstly,  multiresolu¬ 
tion  properties  are  useful  to  quantify  the  relative  contributions  of  the  large  scales  vs.  the 
small  scales.  Secondly,  the  extent  to  which  the  large  scales  dominate  the  aero-optical  in¬ 
teractions  is  useful  for  development  of  methods  for  reduction  of  beam  aberrations  such  as 
the  two  principal  approaches  on  flow  control  summarized  above  and  changes  in  the  flow 
geometry.  In  the  results  presented  below,  we  examine  the  effect  of  modifying  the  resolu¬ 
tion  of  the  refractive-index  held  for  various  resolutions.  These  coarse-grained  fields  are 
then  utilized  to  examine  the  resolution  effect,  i.e.  the  extent  to  which  there  is  sensitivity 
or  robustness,  on  the  aero-optical  interactions  and  aberrations. 

As  a  reminder,  we  note  that  one  of  the  key  features  of  the  work  reported  in  this  grant 
period  is  the  combination  of  an  experimental  approach  with  a  computational  approach  in 
order  to  examine  the  aero-optical  interactions  along  the  beam  propagation  path,  retain¬ 
ing  the  large-Reynolds-number  capability  of  the  experiments  while  including  the  ben¬ 
efits  of  computations.  We  confine  our  discussion  to  the  near-held  aero-optical  distor¬ 
tions  for  which  geometrical  optics  provides  a  useful  approximation  [7],  whereas  far-held 
optical  beam  propagation  requires  Fourier  optics.  The  present  combined  experimen¬ 
tal/computational  approach  can  be  appreciated  by  realizing  that  in  purely-computational 
aero-optics,  [52,  60,  21,  53,  7,  51,  8]  as  in  computational  fluid  dynamics,  there  are  currently 
limitations  regarding  the  Reynolds  numbers  that  are  achievable  even  with  massively- 
parallel  computers  in  direct  numerical  simulation  (DNS)  techniques.  Current  computa¬ 
tions  are  limited  to  moderate  Reynolds  numbers  unless  subgrid-scale  modeling  is  em¬ 
ployed  as  in  large-eddy  simulations  (LES).  In  contrast,  whole-held  imaging  experiments 
such  as  in  the  present  study  at  large  Reynolds  numbers  and  at  Unite  resolution  provide 
large-eddy  experimental  results  free  from  any  subgrid-scale  modeling. 
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Figure  15:  Top  Left:  Experimental  refractive  field  image  at  a  full  resolution  of  ~  1000  x 
1000  pixels  for  the  compressible  separated  shear  layer  with  Re  ~6x  106,  Mc  0.4,  and 
p  ~  3  atm.  This  image  corresponds  to  laser-induced  fluorescence  arising  from  aero-optical 
mixing  effects  as  opposed  to  density  effects.  Top  Right:  Reduction  of  the  resolution  to 
~  75  x  75  pixels  produced  by  coarse  graining  the  original  experimental  high-resolution 
refractive  field  image.  Bottom  Left:  Reduction  of  the  resolution  to  ~  50  x  50  pixels 
produced  by  coarse  graining  the  original  experimental  high-resolution  refractive  field 
image.  Bottom  Right:  Reduction  of  the  resolution  to  ~  30  x  30  pixels  produced  by 
coarse  graining  the  original  experimental  high-resolution  refractive  held  image. 
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Figure  16:  Surface  rendering  of  the  two-dimensional  spatial  refractive  index  field  for  the 
middle  vertical  third  of  the  image  in  the  top  left  of  figure  15,  at  a  spatial  resolution  of 
150  x  600  pixels.  This  surface  shows  the  variations  of  the  refractive  index  across  the 
transverse  instantaneous  extent  of  the  shear  layer  encountered  by  the  laser  wavefronts. 
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Figure  17:  Surface  rendering  of  the  two-dimensional  spatial  refractive  index  field  for  the 
middle  vertical  third  of  the  image  in  the  top  right  of  figure  15,  at  a  spatial  resolution  of 
100  x  400  pixels.  This  surface  shows  the  coarse-grained  variations  of  the  refractive  index 
across  the  transverse  instantaneous  extent  of  the  shear  layer  encountered  by  the  laser 
wavefronts. 
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Figure  18:  Surface  rendering  of  the  two-dimensional  spatial  refractive  index  field  for  the 
middle  vertical  third  of  the  image  in  the  bottom  left  of  figure  15,  at  a  spatial  resolution  of 
50  x  200  pixels.  This  surface  shows  the  coarse-grained  variations  of  the  refractive  index 
across  the  transverse  instantaneous  extent  of  the  shear  layer  encountered  by  the  laser 
wavefronts. 
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Figure  19:  Surface  rendering  of  the  two-dimensional  spatial  refractive  index  field  for  the 
middle  vertical  third  of  the  image  in  the  bottom  right  of  figure  15,  at  a  spatial  resolution 
of  25  x  100  pixels.  This  surface  shows  the  coarse-grained  variations  of  the  refractive 
index  across  the  transverse  instantaneous  extent  of  the  shear  layer  encountered  by  the 
laser  wavefronts. 
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Figure  20:  Surface  rendering  of  the  cumulative  optical  path  difference  variations  for  laser 
wavefront  propagation  through  the  refractive  index  held  of  figure  16  across  the  transverse 
extent  of  the  shear  layer.  This  surface  rendering  is  for  the  full  spatial  resolution  of 
~  1000  x  1000  pixels. 
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Figure  21:  Surface  rendering  of  the  cumulative  optical  path  difference  variations  for  laser 
wavefront  propagation  through  the  coarse-grained  refractive  index  held  of  figure  17  across 
the  transverse  extent  of  the  shear  layer.  This  surface  rendering  is  for  the  reduced  spatial 
resolution  of  ~  75  x  75  pixels. 
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Figure  22:  Surface  rendering  of  the  cumulative  optical  path  difference  variations  for  laser 
wavefront  propagation  through  the  coarse-grained  refractive  index  held  of  figure  18  across 
the  transverse  extent  of  the  shear  layer.  This  surface  rendering  is  for  the  reduced  spatial 
resolution  of  ~  50  x  50  pixels. 
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Figure  23:  Surface  rendering  of  the  cumulative  optical  path  difference  variations  for  laser 
wavefront  propagation  through  the  coarse-grained  refractive  index  held  of  figure  19  across 
the  transverse  extent  of  the  shear  layer.  This  surface  rendering  is  for  the  reduced  spatial 
resolution  of  ~  30  x  30  pixels. 
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Multiresolution  results  for  both  the  refractive-index  behavior  and  the  aero-optical  cu¬ 
mulative  interactions  are  shown  in  the  sequences  in  figure  15  and  figures  16  through  19. 
The  first  sequence,  in  figures  15,  shows  the  coarse-graining  in  refractive  field  information 
as  the  resolution  is  reduced.  The  second  sequence,  in  figures  16  through  19,  shows  vi¬ 
sualizations  of  the  variations  in  the  mixing  signature  of  the  refractive  index  across  the 
shear  layer.  The  third  sequence,  in  figures  20  through  23,  shows  the  multiresolution 
effects  on  the  cumulative  aero-optical  interactions.  It  is  directly  evident,  in  these  three 
sequences,  that  the  primary  structure  of  the  aero-optical  wavefront  aberrations  along  the 
beam  propagation  paths  is  highly  robust  to  the  coarse-graining  or  resolution  reduction. 
This  is  a  key  finding  with  implications  both  experimentally  and  computationally,  vis-a- 
vis  bandwidth  requirements  for  experimental  flow/beam  control  and  resolution  require¬ 
ments  for  computational  simulations  respectively,  regarding  the  aero-optical  interactions 
for  directed-energy  research  as  well  as  development. 

3.  Conclusions  and  Implications 

In  the  research  reported  in  this  grant  period,  a  key  focus  is  the  role  of  refractive  interfaces 
in  large  scale  aero-optical  interactions  in  high  Reynolds  number  compressible  separated 
shear  layers.  Direct  imaging  of  the  refractive  index  field  and  interfaces  is  conducted 
using  laser  induced  fluorescence  of  acetone  vapor  molecularly  premixed  in  air  and  digital 
photography  techniques.  High  resolution  two-dimensional  spatial  images  of  the  refractive 
index  field  have  been  recorded.  The  Reynolds  number  based  on  the  large  scale  visual 
thickness  is  Re  ~  6  x  106,  the  convective  Mach  number  is  Mc  ~  0.4,  and  the  test  section 
pressure  is  p  ~  3  atm.  The  refractive  index  fields  have  been  examined  to  identify  the 
dominant  interfacial  contributions  to  the  propagation  and  distortion  of  optical  wavefronts. 
Highly  irregular  networks  of  spatially-isolated  high  gradient  interfaces  are  evident  both 
in  the  flow  interior  and  at  the  outer  boundaries.  These  high  gradient  interfaces  are  found 
to  be  primarily  responsible  for  the  large  scale  aero-optical  interactions.  Specifically,  a 
physical  model  utilizing  these  high  gradient  boundaries  is  able  to  reproduce  well  the 
large  scale  optical  wavefront  distortions.  Significant  robustness  is  found  for  the  optical 
distortions  in  terms  of  reduction  in  the  refractive  field  information  by  retaining  high 
refractive  gradient  interfaces.  This  has  implications  for  optimization  studies  in  aero- 
optics  because  it  suggests  that  effective  modification  of  the  optical  aberrations  may  be 
sought  in  terms  of  altering  the  structure  of  the  high  gradient  interfaces. 

The  multiresolution  study  based  on  the  combined  experimental/computational  approach 
provides  support  for  ongoing  studies  of  aero-optical  interactions  in  laser  beam  propa¬ 
gation  through  highly-separated  turbulent  flows.  The  robustness  to  resolution  reduc¬ 
tion,  i.e.  coarse-graining,  is  an  encouraging  finding  from  an  applications  point  of  view, 
e.g.  in  terms  of  large-scale  flow  organization  and  disorganization  [61,  14],  The  large 
Reynolds  number  flow  conditions  are  achieved  with  a  specially-designed  experimental 
laboratory  facility.  This  variable-pressure  flow  facility  enables  direct  imaging  of  the  tur¬ 
bulent  refractive-index  field  at  Reynolds  numbers  of  Re  ~  106  in  controlled  laboratory 
flows.  These  large  Reynolds  numbers  are  necessary  in  order  to  exceed,  in  this  case 
by  two  orders  of  magnitude,  the  minimum  value  of  Re  ~  104  representative  of  fully- 
developed  turbulence.  Laser-induced  fluorescence  in  purely-gaseous  flows,  using  acetone 
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vapor  seeded  in  air,  is  utilized  to  directly  image  the  turbulent  refractive-index  spatial 
fluctuations  along  the  beam  propagation  paths  with  a  high-resolution  intensified  CCD 
camera.  Direct  imaging  of  the  refractive-index  fields  is  enabled  by  propagating  pulsed  ul¬ 
traviolet  laser  sheets  through  the  highly-separated  shear  layers.  The  present  methodology 
provides  a  means  to  study  the  laser  beam  wavefronts  along  the  entire  propagation  path  in 
large  Reynolds  number  flows,  which  are  otherwise  inaccessible  for  purely-computational 
studies  that  would  be  limited  to  lower  Reynolds  numbers  or  small-scale  turbulence  mod¬ 
eling.  The  present  multiresolution  study  enables  an  examination  of  the  sensitivity  of 
the  aero-optical  aberrations  with  regard  to  the  structure  of  the  refractive  interfaces  en¬ 
countered  along  the  propagation  path.  The  evidence  of  robustness  of  the  aero-optical 
aberrations  to  the  amount  of  reduction  in  refractive  held  information,  over  a  wide  range, 
provides  a  key  ingredient  for  experimental  and  computational  advances  in  aero-optical 
aspects  of  directed-energy  research  and  development. 
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